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Abstract 

A new CAD program, PolyCAD, designed for X-ray photon tracing in 
polycapillary optics is described. To understand the PolyCAD code and 
its results, the theoretical basis of X-ray transmission by a single cylin- 
drical channel (monocapillary) is discussed first. Then the simplest cases 
of cylindrically and conically shaped polycapillary optics are examined. 
PolyCAD allows any type of X-ray source to be used: an X-ray tube of 
finite beam dimensions or an astrophysical object can be simulated in 
combination with the polycapillary optics. The radiation distribution im- 
ages formed on a screen located at various focal distances are discussed. 
The good agreement of some of the PolyCAD results with those reported 
in earlier papers validate the code. This is the first paper of a series 
dedicated to the development of an exhaustive CAD program, work is in 
progress to develop the code to include other polycapillary-optics shapes, 
such as semi-lenses and full-lenses. 

OCIS: 340.0340, 340.7470, 080.2740, 230.7380, 999.9999. 



1 Introduction 

The use of polycapillary optics to control X-ray beams in analytical X-ray ap- 
paratus for diffraction and fluorescence analyses is becoming increasingly im- 
portant. In the near future polycapillary optics will be widely used in many 
different fields, e.g., aero-space research, medicine, biology, and so on [1] . Our 
strong interest in the development of these devices led to the idea of creating a 
CAD program that allows one to simulate the process of radiation propagation 
through polycapillary optical systems and to visualize the radiation distribu- 
tions at the optical output. At present there are very few procedures available 
for evaluating X-ray transmission by capillary structures in the case of peculiar 
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configurations, sec for instance [2, 3, 4, 5, 6, 7, 8] . One of the first papers [2] 
reported coherent radiation transmission by a hollow glass pipe. This code also 
considered i) the presence of a rough surface as an X-ray anomalous disper- 
sion effect, ii) radiation penetration into the channel wall, and iii) the possible 
presence of micro-dust inside the channel. 

Another algorithm [3] for X-ray transmission by capillaries of various shapes 
runs in the geometrical optics approximation; for this reason is simpler and 
more fiexible than the previous model. 

The first X-ray tracing codes were developed in 1992 [4, 5] and by the In- 
stitute for Roentgen Optics, for a review see [11] . However, these simulations 
used a number of simplifications based either on the cylindrical symmetry for 
a shaped capillary system or on the meridional ray approximation, which is 
valid only for describing radiation propagation through a blend channel or a 
monocapillary concentrator. 

To the authors' knowledge, the most advanced and complete softwares for 
X-ray tracing inside capillary channels to date are three: the first code [6, 9] 
traces the trajectory of each photon including the corrections of absorption 
and roughness; in according with some experimental results, the authors could 
obtain the roughness of the channel. The second one [7] uses the SHADOW ray 
tracing software, adapting it to channel shape. 

The last code [8, 10] uses a Monte Carlo simulation for X-ray radiation prop- 
agation through hollow channels. The theoretical results obtained by this code 
agree quite well with the experimental data, although the algorithm is rather 
simple due to its geometrical optics approximation and to the circular cross 
section of the channel shape. 

The main aim of practically all the simulation codes is to optimize both the 
channel size and the optical shape in order to obtain highly efficient optical sys- 
tems. Obviously, this is very important from the viewpoint of the development 
of capillary optics technology. However, analysis of the radiation distribution 
features, reported in a number of papers [12, 13, 14] is particularly interesting. 
(Some of the earlier publications are cited in reference [14]). In these papers 
the experimental data were validated by means of analytical estimations based 
on the wave theory of radiation propagation. 

Before we continue any further, let us clarify the terminology used in this 
work: 1) a monocapillary is a single capillary; 2) a polycapillary is a set of 
closely and packed monocapillaries; 3) a lens is a device that concentrate the 
radiation in a point or a small region; 4) a semi-lens, or a half-lens is a device 
that can concentrate a quasi-parallel beam in a point and vice versa in the 
reverse geometry. 

Here, we introduce a new X-ray tracing code for polycapillary optics, named 
PolyCAD. A previous version was designed for cyhndrical optics only [15, 16] . 
Now, the software can simulate monocapillary and polycapillary optics with 
cylindrical and conical shapes, while work is in progress to include lenses and 
semi-lenses. The advantage of the code lies in its precise mathematical solutions 
for each given optical shape. Comparison of the results of PolyCAD and of the 
previous algorithms [5, 6, 7, 8, 11] revealed some differences, due to the fact that 
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PolyCAD is free of many of the algorithmic constraints. 

In the first part of our paper, the theoretical basis of a cylindrical monocap- 
illary is explained; in the second, the numerical results and the simulations are 
reported both for cylindrical and for conical mono- and polycapillaries, using a 
point-source and an extended-sources, like a conventional X-ray tube. 



2 Theoretical ground 

Thanks to some simple formulas it is possible to estimate the correct behavior 
of the rays inside the capillary. The total intensity 1(6) , evaluated at the angle 
9 behind a capillary, is given by all the contributions of the photon families 
that pass through the capillary; each family is defined by the number of its 
reflections: 

1(6) = io{9) + h{e) + --- + iNie) (1) 

where N is the maximum number of reflections. 

It is clear that the intensity lm{0) is strictly dependent on the reflection 
coefficient Rq, namely, on the m*'' power of Rq: 

loie) = Roe€i-eo,Oo) (2a) 



In{0) k R^9e{-9N,eN) (2b) 

This contribution is drastically small if the number of reflections is high, 
because Rq is always less than one, (Example - Rq = {0.99, 0.8, 0.7} — > i?Q° = 
{0.61, 1.42 X 10-^ 1.80 X 10-8}). 

In order to consider the total radiation intensity on a screen behind a mono- 
capillary, I{9) has to be integrated over the whole plane, so: 

itotie) = J J d^^dd'^ (X /^<^|e ^^(^)} = 2^ |e (3) 

where the last mathematical passage is possible due to the axial symmetry of 
the system. 

As can be seen from Fig. 1, if the screen is placed at PI, near the focal point 
/2, characterized by the ray family with a maximum number of two reflections, 
the image on the screen is smaller than in any other configuration. Moreover, 
for only two reflections, there are flve different areas of interaction (Figs. 1 and 
4). However, if the screen is placed at P2, far from the focal point, the picture 
on the screen does not have a simple shape, but has some zones that are more 
intense than others (Figs. 1 and 6). 

Now that the radiation intensity has been defined (Eq. 3), the inner surface 
of the capillary can be divided into a discrete set of zones according to the 
number (1, 2,. . . ) of refiections of each photon trajectory. If the origin of the 
reference system {x, y, z) is located at the end of the capillary and the z-axis 
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coincides with the capillary axis, each zone along the surface has a length given 
by the following equation: 



z{e) ^ fx i{0), {0 « 1) 



(4) 



where / is the distance between the end of the capillary and the screen. Looking 
at Fig. 2 and considering only the upper propagation case, each area has a 
maximum angle 9 and a portion of interaction on the screen, defined by: 
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where d is the diameter of the capillary and i?™ is the m*'* reflection coefficient. 
The maximum angle for each reflection family can be evaluated as: 
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From these formulas it is clear that the difference between two consecutive 
maximum angles is a constant and depends only on the properties of the system. 
When these properties are flxed, the maximum number of reflections is simple 
to evaluate: 
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Since there must be at least one reflection, Eq. 7 becomes the well known 
expression: 

L " 
2s 



N = 



+ 1 



(8) 
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Equation 8 states that the maximum number of reflections depends only 
on the length and on the diameter of the capillary; this result is exact if the 
source stays along the monocapillary z-axis. Equation 8 also remains a good 
approximation if the source is close to the optical axis (Fig. 12). Obviously, 
in the general case, there are some problems because an off-axis source makes 
some zones inside the capillary unsuitable for photon reflection [17] . 



3 The code 

In our previous paper [15] the code was based only on the geometrical symmetry 
of the system, so the radiation image on a screen could be evaluated only if the 
shape of capillary was cylindrical. This new version of PolyCAD is completely 
different from the old one. 

While upgrading the code to include other capillary shapes, the complicated 
definition of the reflection angle a^e/ introduced some problems. For a cylin- 
drical capillary it is easy to find the following relation for are/: 

/cosQ!\ f COS a 1 , . 

aref = arccos = arccos < f ; : — --^ } (9) 

V cos a; / |^ cos [arctan (tan a x sin t*)] J 

where a is the zenithal angle and 6 is the angle on the inlet plane (see Fig. 3 
and Ref. [15]). 

In the conical case, the u> angle, formed by the two projections on the tangent 
plane of the X-ray photon trajectory and of the z-axis, is still a function of a 
and 9, but now there is also an indirect dependence on the conical semi-opening 
angle /?. This indirect dependence is due to the fact that the trajectory of a 
photon is generally oblique with respect to the capillary axis. Thus, a new 
parameter that takes into account the off-axis angle must be included, which 
does not allow us to transfer in a simple way a 3-D problem into a 2-D one, as 
done in the case of cylindrical optics [15] . Nevertheless, by knowing the photon 
starting point, i.e., the source point Pg and the inlet point Pq, it is possible 
to define the vector-direction of the photon. An equation system between the 
photon path and the capillary surface equation makes it possible to describe the 
trajectory of the photons inside any generic capillary optics. 



4 Numerical results 

To show how PolyCAD works, we report the most significant results. We would 
like to emphasize that PolyCAD can simulate any kind of source, from a point 
source, located at finite or infinite distance, to an extended source, and the 
latter can also have a 3-D shape. 

4.1 Cylindrical monocapillary 

We first considered a monocapillary or polycapillary with cylindrical channels. 
Even though the system geometry in this case is so simple that the results can 
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be calculated manually, to check that the software worked well the results were 
compared with those of the old version of PolyCAD designed just for cylindrical 
capillaries [15, 16] . 

As a preliminary step we considered a point source located at finite distance. 
The amplitude of the incident angle a can be chosen randomly from zero to any 
prefixed value. Obviously, in the best case this angle is equal to the critical 
angle, a = O,.. Moreover, each X-ray photon will have an inlet angle S in the Ip 
plane that is strictly connected with the intersection point Pq {xq, yo) (Fig. 3). 

To simplify the analysis of the radiation distribution behind the capillary 
system, the image shape was considered in two different screen positions: i) one 
placed at the focus f2 of the monocapillary or polycapillary and ii) the other 
beyond the f2 position (Fig. 1). 

To compare these calculations with the results obtained by Dabagov and 
Marcelli [12] , the length of the cylindrical monocapillary was chosen such as 
to provide at least a double reflection mode of propagation. The final results 
obtained are illustrated in Figs. 4, 6 and 8. 

The conditions for the simulation are: i) a 1 keV point source and ii) a 
cylindrical monocapillary 10 cm long with a radius of p = 10~^ cm. In Figs. 4 
and 5 a single central spot is present and four rings formed by the convergent 
rays related to the different reflection numbers (compare Figs. 4 and 1). 

In Figs. 6 and 7 many sharp rings are present. The central ring is given 
by the X-ray photons that cross the monocapillary without interacting with the 
surface. The second is formed by X-ray photons that hit, once at worst, the 
monocapillary inner surface. The third ring is due to X-rays that produce only 
one reflection with the channel. The fourth ring is due to X-rays that make one 
or two reflections. Finally the outer ring is formed by ray vectors that undergo 
only two reflections. 

4.2 Cylindrical polycapillary 

In the case of a cylindrical polycapillary the simulation parameters are i) a 1 
keV point source; ii) a polycapillary 10 cm long with a radius of p = 1 cm; iii) 
a single-channel radius of 10^^ cm. 

Looking at Fig. 8, for this configuration it is possible to observe a slightly 
bigger spot in the center, while the point density distribution of the halo de- 
creases going from the center to the periphery. This last effect is due to the fact 
that each X-ray photon interacts many times with the channel surface and the 
number of interactions increases going from the center to the periphery. 

4.3 Conical monocapillary 

As the next step we will consider a conical capillary with an X-ray point source 
located at finite distance. Even though the amplitude of the incident angle 
a can be chosen in a random way from zero to any prefixed value, in order to 
maximize the number of incoming photons that undergo total external refiection 
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the following relation should be satisfied: a+/3 = 6c, where /3 is the semi-opening 
angle of the cone. 

It is not a simple job to define the focal plane for a conical capillary be- 
cause each photon after each interaction with the channel surface undergoes a 
number of refiections with increasing incidence angular values: a + 2(3, a + A(5, 
. . .etc., according to the different directions of the incident rays. Taking this 
into account wc decided to choose the position of the focal plane by minimizing 
the divergence of the exit beam. To compare the results of the conical and 
cylindrical capillaries, we again chose X-ray photons that undergo at least two 
reflections. 

The parameters of the conical monocapillary are similar to the cylindrical, 
i.e., i) the capillary length is 10 cm; ii) the entrance radius pi = 10~^ cm, while 
the exit radius is p/ = 0.8 x 10^"'^ cm; iii) the energy of the photon is still 1 keV. 

As for the cylindrical capihary, simulations of the focal and out-of-focus 
radiation distributions are shown in Figs. 9, 10 for the conical monocapillary. 

4.4 Conical polycapillary 

For a conical polycapillary the simulation results are shown in Fig. 11. The 
parameters are: i) a 1 keV point source; ii) a polycapillary 10 cm long with 
the entering radius of pi = 1 cm, while the exit radius is pf = 0.8 cm; iii) a 
single-channel radius of 10~^ cm and 0.8 x 10~^ cm, respectively for the entering 
and the exit plane. 

In a following paper we will discuss in more detail the differences between 
conical and cylindrical capillaries. In particular we will show that the conical 
capillary has a more intense central spot and also a bigger halo. 

4.5 Cylindrical monocapillary, off axis case 

In this section, we discuss the problem of an X-ray source located off-axis. The 
point source can be located at finite or at infinite distance from the capihary, 
so the behavior of a cylindrical monocapillary system differs accordingly (Figs. 
12 and 13). 

Figure 12 shows a simulation of a cylindrical monocapillary with a point 
source that is not placed along the z-axis, with coordinates: Pg = (0.05, 0.05, 3.3) 
cm. The other physical parameters are capillary length = 10 cm, capillary radius 
= 10~^ cm. 

We would like to point out that in this configuration the image has a sym- 
metric shape with respect to an ideal line, which has a 45° angular inclination 
due to the mirror-like behavior of each half of the channel. A large circular 
halo is also present, due to all the photons that cross the capillary without any 
interaction. 

In Fig. 13 the situation is completely different because the point source is 
again located off-axis but at an infinite distance from the capillary. The zenithal 
and equatorial angles of the beam are respectively {a = 1.3°, 5 = 45°), while 
the channel parameters are, as usual, length = 10 cm, radius = 10~^ cm. This 
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configuration means that some parts of the internal channel surface cannot be 
reached by X-ray photons, so now certain regions in the image are completely 
empty. 

It is essential to emphasize that the zenithal angle is so large that many 
regions of the capillary wall will not be active for X-ray propagation. 

4.6 Extended sources 

The power of PolyCAD program is that it can deal with any geometrical and 
optical configuration. It means that we may easily simulate the behavior of a 
source of finite dimensions like, for example, conventional X-ray tube. In order 
to treat this situation an offset for a single point-source is introduced in the 
program. In such a mode a finite area source can be simulated by random 
distribution of points inside a specific area. As practical application let to 
consider a 'Long Fine Focus' of X-ray tube (for instance, a Cu anode, 8keV), 
where the electron spot on the anticathode has the dimension of 0.4 x 12 mm,. 
Such kind of tube allows us to have both an 'optical point-focus' (Fig. 14) and 
an 'optical hne-focus' (Fig. 15), where in these simulations we used cylindrical 
polycapillary optics. The source dimensions are respectively 0.4 x 1.2 mm and 
0.04 X 12 mm. 

5 Conclusion 

In this paper we have reported the first part of a PolyCAD program designed 
for capillary optics. Using the ray optics approximation in the ideal case of total 
external refiection, i.e., without absorption effects by the optical channel walls, 
this algorithm allows us to i) simulate the passage of the X-ray beam inside the 
capihary channels for sources of different shapes, ii) visuahze the spot images 
formed on a screen at different distances from the optics. 

After a description of the theoretical basis and of the computational details, 
we reported the numerical analysis for various source-optics configurations. In 
the case of a cylindrical monocapillary, we found good agreement with previous 
results [12] , both for single and for multiple reflection modes. 

We would hke to point out that this code can accurately describe the radia- 
tion distribution behind the optics, as shown in Figs. 5 and 7. When the source 
is located along the axis of the capillary, and in particular taking into account 
a single channel, photons are present in all the image areas; however it is clear 
that by changing the source-optics conditions, the photons could be absent in 
some areas. 

At present we are working on increasing the PolyCAD program capabilities 
to deal with other polycapillary-optics shapes, such as semi-lenses and full- 
lenses. However, to understand the fine features of X-ray propagation through 
polycapillary lenses it is mandatory to consider the X-ray wave interaction with 
the inner capihary surface (see details in the review [14]). This will constitute 
the future development of PolyCAD. 
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List of Figure Captions 



Fig. 1. Cross section of cylindrical monocapillary with on-axis source. The 
expected beam intensity on a screen located in two different positions is shown. 
The image shape changes according to position, but note that areas without any 
counts are not allowed. Focal points /I and /2 are due to rays that undergo 
one and two reflections, respectively. 

Fig. 2. Length of capillary L\ distance from source s; distance from screen /; 
channel radius d. The maximum angles allowed for each ray family are defined. 
Rays entering between zq and zi undergo only one reflection, rays between zi 
and Z2 undergo only two, and so on. 

Fig. 3. Three-dimensional drawing of photon reflection from a cylindrical and 
a conical capillary wall. Bold line shows photon path. Ip is the transverse cross 
section plane at the capillary entrance. See text for symbol meaning. 
Fig. 4. Cylindrical monocapillary. Image of a point source on a screen located 
at /2 (Cf. Fig. 1. The intense central spot is formed by all the ray families. 
Channel dimensions are radius = 0.1 cm and length L = 10 cm. As per theory 
there are no empty areas in the image. 

Fig. 5. Three-dimensional density distribution of the photons in Fig. 4. Source 
and screen are in focal points s and /2, respectively. The radiation distribution 
inside the spot is never equal to zero. 

Fig. 6. Point-source image on a screen located at P2 (Cf. Fig. 1) out of focus 
behind the cylindrical channel. The distribution consists of concentric rings of 
different intensities with no empty areas. 

Fig. 7. Three-dimensional density distribution of the photons in Fig. 6. Source 
and screen are in points s and P2, respectively. A section taken along the y-axis 
conflrms that even in the center of the spot the radiation distribution is never 
equal to zero. 

Fig. 8. Image of a point source on a screen behind a cylindrical poly capillary. 
The main difference between the polycapillary and monocapillary is that the 
halo is not homogeneous, Polycapillary parameters: length 10 cm and radius 1 
cm. Radius of each capillary is 0.001 cm. 

Fig. 9. Image formed by a conical monocapillary in the best geometrical condi- 
tions with the screen on the focal point (see text for details). Channel parame- 
ters: radius at capillary entrance 10~^ cm; radius at exit 0.8 x 10~^ cm; energy 
1 keV; monocapillary length L=10 cm. 

Fig. 10. Conical monocapillary output image on a screen located out of focus. 
Here, the concentric rings have higher intensity than the central spot. Parame- 
ters same as in Fig. 9. 

Fig. 11. Point-source image on a screen behind a conical polycapillary. As in 
Fig. 8 the main difference is that the halo is not homogeneous. Polycapillary 
parameters: length = 10 cm, entrance radius 1 cm and exit radius 0.8 cm. 
Radius of each capillary is 10~^ cm at entrance and 0.8 x 10~^ cm at exit. 
Fig. 12. Image due to a single cylindrical channel, when point source off-axis. 
Note that there are no areas without any count. Capillary parameters same as 
in Fig. 4. Source is at point (0.05, 0.05, 3, 3) cm from entrance plane. 
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Fig. 13. Image with point source at infinite distance (a typical astrophysics 
source). Cylindrical monocapillary parameters same as in Fig. 12. Source 
parallel beam has zenithal angle a = 1.3° and equatorial angle 6 = 45° (Cf. 
Fig. 3). 

Fig. 14. The image formed on the screen by the cylindrical polycapillary, when 
the 'optical point focus' of a Cu tube has been used. 

Fig. 15. The image behind cylindrical polycapillary sample when the 'optical 
fine focus' of a Cu tube has been used. 
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